Nonlinear adiabatic passage from fermion atoms to boson molecules 



E. Pazy 1 , I Tikhonenkov 1 , Y. B. Band 1 , M. Fleischhaucr 2 , and A. Vardi 1 
1 Department of Chemistry, Ben-Gurion University of the Negev, P.O.B. 653, Beer-Sheva 84105, Israel and 
Fachbereich Physik, Technische Universitat Kaiserslautern, D67663, Kaiserslautern, Germany 

We study the dynamics of an adiabatic sweep through a Feshbach resonance in a quantum gas 
of fermionic atoms. Analysis of the dynamical equations, supported by mean-field and many-body 
numerical results, shows that the dependence of the remaining atomic fraction T on the sweep 
rate a varies from exponential Landau-Zener behavior for a single pair of particles to a power-law 
dependence for large particle number N. The power-law is linear, Faa, when the initial molecular 
fraction is smaller than the 1/N quantum fluctuations, and F oc a 1 ^ 3 when it is larger. Experimental 
data agree better with a linear dependence than with an exponential Landau-Zener fit, indicating 
that many-body effects are significant in the atom-molecule conversion process. 



Adiabatic evolution is an important tool for quantum 
state engineering. The adiabatic theorem ensures that 
an initial nondegenerate eigenstate remains an instanta- 
neous eigenstate when the Hamiltonian changes slowly. 
When eigenstates become nearly degenerate, the Landau- 
Zener (LZ) model Q is a paradigm for explaining how 
transitions occur. 

Adiabatic sweeps across an atom-molecule Feshbach 
resonance have recently been used to convert degenerate 
fermionic atomic gases containing two different internal 
spin states to bosonic dimer molecules Forma- 
tion of a molecular condensate has also been observed us- 
ing both adiabatic sweeps and three-body recombination 
processes Q. In this Letter we show that for adiabatic 
Feshbach sweeps that convert degenerate fermionic atoms 
to diatomic molecules, the LZ behavior for a single pair 
of particles is dramatically changed due to many-body ef- 
fects. The fraction of unconverted atoms is shown to be 
a power-law in the sweep rate, rather than exponentially 
small as predicted by an essentially single-particle, linear 
LZ model Q, Q • The exact power-law is determined by 
the significance of quantum fluctuations. En route to this 
result we also find that, for a ladder of atomic states filled 
by fermionic atoms, the atom-molecule sweep efficiency is 
unaffected by atomic dispersion, and all fermionic atoms 
can go over to molecules, in contrast to the linear LZ 
model. 

We consider the collisionless, interaction representa- 
tion, single bosonic mode Hamiltonian 0, 0, 0, 0, 0, 

H = }] ekC^ „Ck,o- + £(t)6p6o 



+9 52 c k>T c_ kj ;6j + B.c. 



(1) 



where e k = h 2 k 2 /2m is the kinetic energy of an atom with 
mass m, and g is the atom-molecule coupling strength. 
The molecular energy £(t) = at is linearly swept at a rate 
a through resonance to induce adiabatic conversion of 
Fermi atoms to Bose molecules. The annihilation opera- 
tors for the atoms, Ck.tr, obey fermionic anticommutation 
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FIG. 1: Many-body collective dynamics of adiabatic passage 
from a fermionic atomic gas into a molecular BEC for five 
pairs of fermionic atoms, (a) Sweep rate a — 2g 2 N, (b) Sweep 
rate a — f/ 2 JV/4. Overall efficiency is independent of atomic 
dispersion in both (a) and (b). 



relations, whereas the molecular annihilation operator bo 
obeys a bosonic commutation relation. 

We find that, provided that all atomic levels are swept 
through, the adiabatic conversion efficiency is completely 
insensitive to the details of the atomic dispersion. Fig. 
shows exact numerical results for the adiabatic conver- 
sion of five atom pairs into molecules, for different values 
of the atomic level spacing (and hence of the Fermi en- 
ergy Ep). It is evident that, while the exact dynamics 
depends on Ep, levels are sequentially crossed, leading 
to the same final efficiency regardless of the atomic mo- 
tional timescale. In particular, in the limit as a — > 
it is possible to convert all atom pairs into molecules. 
This is a unique feature of the nonlinear parametric cou- 
pling between atoms and molecules, which should be con- 
trasted with a marginal conversion efficiency expected for 
linear coupling. Since the exact energies e k do not af- 
fect the final fraction of molecules, we use a degenerate 
mode l FH| 111 ITij with e k = e for all k. In the spirit of 
Refs. we define the operators: 
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Ck cr is the conserved to- 



wherc N 

tal number of particles. It is important to note that 
J-,J+ 7 J Z do not span SU{2) as [J + ,JJ\ is a quadratic 
polynomial in J z . We also define J x = J + + J_ and 
J y = —i(J+ — J-). Up to a c- number term, Hamilto- 
nian (JTJ takes the form 



N I 

H=- I A(t)J z 
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where A(t) = 2e — £(t). Defining a rcscalcd time r = 
VA" gt, we obtain the Hcisenberg equations of motion for 
the association of a quantum-degenerate gas of fermions, 
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which depend on the single parameter 8(t) = 
A(t)/y/Ng = (a/g 2 N)r. We note parenthetically that 
precisely the same set of equations, with J z — > — J z and 
q — > p/2, is obtained for a two-mode atom-molecule BEC 
|l5| , providing another perspective on the recently ob- 
served mapping between the two systems [lj, [lj, [l4| . 

We first consider the mean-field limit of Eqs. (QJ, re- 
placing J x ,J y , and J z by their expectation values it, 
i;, and io which correspond to the real and imaginary 
parts of the atom-molecule coherence and the atom- 
molecule population imbalance, respectively, and omit- 
ting the quantum noise term y/2(l + J z )/N. In this limit, 
the equations depict the motion of a generalized Bloch 
vector on a two-dimensional surface, determined by the 
conservation law, 



u 2 +v 2 = ^(w- 1) 2 (> + 1). 



(5) 



Hamiltonian J2J is then replaced by the classical form 



gN^ 2 



(Sw + y/(l + w)(l -w 2 )cos6^ 



H(w,6;A) 

'(6) 

with 9 = arctan(w/u). 

To study the atom-molecule adiabatic passage, we 
closely follow the method of Ref. ^(|. The eigenvalues 
of the atom-molecule system at any given value of 5 cor- 
respond to the fixed points (uq,vq,wo) of the classical 
Hamiltonian JJJJ or the mean-field limit of Eqs. 
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(w - 1) (3w Q + 1) = 8u Q . 
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FIG. 2: Equal-energy contours of Hamiltonian @ plotted as 
a function of w and 6 for different detunings S. w = 1 is all 
atoms and w = — 1 is all molecules. The various fixed points 
corresponding to adiabatic eigenvectors are marked by (blue) 
squares, (red) circles and (black) stars. 



The number of fixed points depends on the parameter S. 
The point uq = vq = 0, wq = 1 is stationary for any value 
of S. Using Eqs. |J5J and 10, other fixed points satisfy 



(3ui 



i(wo + 1) 



= S 2 
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In Fig. |2]we plot phase-space trajectories, corresponding 
to equal-energy contours of Hamiltonian @, for differ- 
ent values of 5. As expected from the plots have the 
symmetry (w, 0; S) <-> (w, 9+n] —5). For sufficiently large 
detuning, \S\ > y/2, Eq. © has only one solution in the 
range — 1 < Wo < 1. Therefore, there are only two (ellip- 
tic) fixed points, denoted by a red circle corresponding 
to the solution of Eq. ©, and a blue square at (0,0,1). 
As the detuning is changed, one of these fixed points 
(red circle) smoothly moves from all-molecules towards 
the atomic mode. At detuning 5 = —\J~2~ a homoclinic 
orbit appears through the point (0, 0, 1) which bifurcates 
into an unstable (hyperbolic) fixed point (black star) re- 
maining on the atomic mode, and an elliptic fixed point 
(blue square) which starts moving towards the molecular 
mode. Consequently, in the regime |5| < y/2 there are 
two elliptic fixed points and one hyperbolic fixed point, 
corresponding to the unstable all-atoms mode. Another 
crossing occurs at 5 = v2 when the fixed point which 
started near the molecular mode (red circle) coalesces 
with the all-atoms mode (black star). 
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The frequency of small periodic orbits around the fixed 
points, f^o, is found by linearization of the dynamical 
equations Q about (u ,vq,wq) and using |JSJ to obtain 



Q 



= ^ 2 + (l-3^o) = J (1 '7 )(3 "° + 5) 
V 4(io + l) 



(9) 



Hence, for \5\ < \[2 the period of the homoclinic trajec- 
tory beginning at (0, 0, 1) diverges. 

Transforming w, 8 into action-angle variables /, 4>, the 
non-adiabatic tunneling probability T at any finite sweep 
rate a is related to the action / accumulated during the 
sweep HE E3, 
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(10) 



where R(I , <fi) is related to the generating function of the 
canonical transformation w, 9 — > I, <t>. We note that, un- 
like the linear pj or Josephson [la. \v\ case, where the 
tunneling probability is linearly proportional to the ac- 
tion increment AI, our choice of variables © causes the 
atomic population at the end of the sweep (and hence, 
r) to be proportional to the square root of AI (since 



\t f )+V 2 (t f )oC 



, where n^ a {tf) is the 



Ek, CT n k, CT (*/) 

population in state |k, a) at the final time tf). Equa- 
tion 1)1 0(1 depicts the familiar rule that in order to at- 
tain adiabaticity, the rate of change of the adiabatic fixed 
points through the variation of the adiabatic parameter 
A, R(I, <f>) A, should be slow with respect to the charac- 
teristic precession frequency <j) = Qq about these station- 
ary vectors. For an adiabatic process where A/<fi — > 0, 
the action (which is proportional to the surface-area en- 
closed within the periodic orbit) is an adiabatic invariant, 
so a zero-action elliptic fixed point evolves into a simi- 
lar point trajectory. Action is accumulated mainly in the 
vicinity of singularities where 4> — f^o —> 0. For linear adi- 
abatic passage Q , such singular points lie exclusively off 
the real axis, leading to exponential LZ transition prob- 
abilities. However, when nonlinearities are dominant, as 
in the Mott-insulating Josephson case 0, 0] and our 
case, there are real singularities, leading to power-law 
dependence of the transfer efficiency on the sweep rate. 

It is clear from Eq. @ that, for the atom-molecule 
conversion problem, a real singularity of the integrand in 
Eq. 1(10(1 exists at wq = 1, where the frequency vanishes 
as Slo ~ g^N(l — wq)- Thus, most the the nonadiabatic 
correction is accumulated in the vicinity of this point (all- 
atoms for fcrmions and all- molecules for bosons). Taking 
the derivative of Eq. JBJl with respect to time, we find 
that the response of the fixed-point velocity to a linear 
sweep rate a is, 
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FIG. 3: Fraction of remnant atoms versus inverse ramp speed 
1/B across the 543 G resonance in a two-component Fermi 
gas of 6 Li. The experimental data (red squares) of Ref. Q, 
which saturates at a remnant of 1/2 [20I . and the mean-field 
calculations obey a linear dependence on sweep rate beyond 
0.4 ms/G. is multiplied by 0.4 ms/G to scale the abscissa 
for the calculated results. The insert shows the best exponen- 
tial fit of the data as a dashed line. 



Having found wo, we can now find the action- angle vari- 
able 4> hi terms of wq: <f> = j <j)dt = j Qo^^-. In the 

vicinity of the singularity we have £Iq ~ g^J N(\ — wq) 
and wq ~ y/2a/gV~N, resulting in 



a 6 



(12) 



Using Eq. 1(12(1 . we finally find that near the singularity, 
(f> — £Iq w g^J N(l — wq) is given in terms of (j) as 



(13) 



Substituting ljT3|) and A — a into Eq. ljTU)) we find that 
the nonadiabatic correction depends on a as 




(14) 



So far, we have neglected the effect of quantum fluc- 
tuations, which are partially accounted for by the source 
term (v / 2/A r )(l + J z ) in Eqs. (gj. As a result, we found 
that Wq does not vanish as wq approaches 1. Conse- 
quently, the remaining atomic population is expected to 
scale as the cubic root of the sweep rate if the initial 
average molecular fraction is larger than the quantum 
noise. However, starting purely with fcrmion atoms (or 
with molecules made of bosonic atoms) , fluctuations will 
initially dominate the conversion process. Equation (0 
should then be replaced by 
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FIG. 4: Many-body calculations for the fraction of remnant 
atoms versus dimensionless inverse sweep rate for various par- 
ticle numbers in the range N = 2 to 800. The many body re- 
sults for large number of particles converge to the mean-field 
results (solid green line) of Fig. [3] 



demonstrating that our previous treatment around wo = 
1 is only valid provided that \wo(U) — 1| 3> l/N. For 
smaller initial molecular population, Eq. should be 
replaced by 



Wq 
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4(w + l) 3/2 N(w + ly^iwo - l) 2 

(16) 

Hence, in the vicinity of wo = 1 the eigenvec- 
tor velocity in the w direction vanishes as wq = 
(y/Na/ gy/8) (wq — l) 2 . The characteristic frequency </> 
is now proportional to (a</>) _1 instead of Eq. I|13|) so that 
A/oca 2 , and G3D1 



r oc a . 



(17) 



Equations l|17|) and (|14fl constitute the main results of 
this work. We predict that the remnant atomic fraction 
in adiabatic Feshbach sweep experiments will scale as 
a power-law with sweep rate due to the curve crossing 
in the nonlinear case. The dependence is expected to 
be linear if the initial molecular population is below the 
quantum- noise level (i.e., when 1 — lUo(ij) *C 1/^0 j an< i 
cubic-root when fluctuations can be neglected (i.e. for 1 — 
wo{U) 3> We note that a similar linear dependence 

was predicted for adiabatic passage from bosonic atoms 
into a molecular BEC 01 ■ 

The analytical predictions illustrated above are con- 
firmed by numerical simulations and existing experimen- 
tal data. Fig. [21 shows the fraction of remnant atoms 

2 

r versus inverse sweep rate computed numerically 
from the mean- field limit of Eqs. (@J. The log-log plot 
highlights the power law dependence obtained beyond 



0.4 ms/G (i.e., for > 1). The numerical results com- 
pare quite well with the experimental data of Ref. Q (red 
squares in Fig- EJ) , indicating that nonlinear effects do in- 
deed dominate the adiabatic passage dynamics in these 
experiments. The power-law fit of the experimental re- 
sults is contrasted with an exponential LZ fit (insert of 
Fig. |3J) which fails to provide an accurate description of 
the observed dependence of efficiency on sweep rate. 

To go beyond the mean-field analysis we carry out ex- 
act many-body calculations of efficiency vs. dimension- 

2 

less inverse sweep-rate, using the methodology of 
[T^ |. The results are shown in Fig. 0] For a single pair 
of particles, N = 2, the quantum association problem is 
formally identical to the linear LZ paradigm, leading to 
an exponential dependence of the remnant atomic frac- 
tion on sweep rate (see insert of Fig.0J). However, as the 
number of particles increases, many-body effects come 
into play, and there is a smooth transition to a power- 
law behavior in the regime a < g 2 N . We note that this 
is precisely the regime where Eq. <|10[) can be used to 
estimate Al and T Q]. 

In summary, we have shown that nonlinear many-body 
effects play a significant role in the atom-molecule con- 
version process for degenerate fcrmionic atomic gases, 
modifying the LZ exponential dependence on sweep rate. 
Experimental data seem to back up this conclusion. 
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